Drifting diffusion on a circle as continuous limit of a multiurn Ehrenfest model 
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which can be rewritten as 
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Tni} s /N, N is the total number of the 
, s denotes the number of balls in the zth 



where fi(s) 
balls, and (m 

urn after s steps. Adding [/;+i(s — 1) — fi-i(s — 1)]/2JV 
to both sides of Eq. (J2J) , we get 
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We study the continuous limit of a multibox Erhenfest urn model proposed before by the authors. 
The evolution of the resulting continuous system is governed by a differential equation, which de- 
scribes a diffusion process on a circle with a nonzero drifting velocity. The short time behavior of 
this diffusion process is obtained directly by solving the equation, while the long time behavior is 
derived using the Poisson summation formula. They reproduce the previous results in the large M 
(number of boxes) limit. We also discuss the connection between this diffusion equation and the 
Schrodinger equation of some quantum mechanical problems. 

PACS numbers: 05.30.-d 

In a previous study we proposed a generalized 
Ehrenfest urn model Q| with N balls and M urns that ar- 
ranged periodically along a circle. The evolution of the 
system is governed by a directed stochastic operation. 
Using the standard matrix diagonalization procedures to- 
gether with a multi-variable generating function method, 
we have solved the problem completely. We found that 
for a generic M > 2 case the average number of balls 
in a certain urn oscillates several times before it reaches 
a stationary value. We also obtained the Poincare cycle 
0, i.e., the average time interval required for the system 
to return to its initial configuration. The result is sim- 
ply given by M N , which indicates that the fundamental 
assumption of statistical mechanics holds in this system. 
Taking M — 2, our model reproduces all the results of 
the original Erhenfest urn model 0. 

In this paper, we further study the continuous limit 
(the large M and TV limit) of the proposed multiurn 
model. We show that by defining a density function p 
as the continuous limit of the fraction f t = (rn^/N, i.e., 
the average number of balls in the ith urn devided by 
N, the continuous limit of the model exists if we also 
define the drifting velocity and diffusion constant appro- 
priately. The evolution of p in spacetime is then governed 
by a differential equation, which can be solved under 
proper initial condition and boundary conditions. The 
results obtained in this paper are in agreement with those 
obtained before by the standard matrix diagonalization 
method. Since for even a generic M-urn and ./V-ball case 
the Poincar'e cycle M is too huge to be experienced, 
the evolution of the system can in practice be treated as 
unrepeatable, thus the average quantities considered here 
become more important than those of microstate details. 
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where At represents the time interval in one step, and Ax 
stands for the center-center distance between two neigh- 
boring urns. Taking the continuous limit, we obtain 
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dt 
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where we have used the substitutions: 
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It is clear that Eq. Q is a diffusion equation. Since the 
model is defined on a circle, we replace x by <f>, v by u>, 
Ax by 8, and the diffusion equation becomes 



dp 
dt 
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We start from the Eq. (4) of Ref.Q: 
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Before further exploring Eq.©, here we give a simple 
and general derivation of the diffusion equation. Note 
that the conservation of probability implies 



dp 
dt 



V - J, 
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where p{r,t) is the probability density and J(r, t) is the 
probability current density. Now, the probability current 



2 



can be written as the sum of two terms, one for the "dif- 
fusion part" , and the other for the "drifting part" of the 
probability carriers (the balls). That is 



M = 30 



-DVp + pv, 



(8) 



where D is the diffusion constant and v is the drifting 
velocity caused by some pumping force. 
Substitute JHJ) into JJJ, we obtain 



(9) 



We further assume that V • v = (incompressible fluid; 
one special case is that v = constant), then we have 



^=DV 2 p-vVp, 



(10) 



which is the desired diffusion equation and has the same 
form as Eq. Q and ©■ 

On a straight line, the above equation becomes Eq. 
and we adopt the boundary condition 



p(oo, t) = p(— oo, t) = 0. 



(11) 



On a circle, Eq. (|10|> becomes Eq. ©, with boundary 
condition 



p(<M)=/# + 27T,i). 



(12) 




M = 60 




Now we find the solutions p for the ID diffusion equa- 
tions on a straight line (QJ and on a circle @), respec- 
tively. Assuming the initially condition: 



FIG. 1: Fraction curves calculated from the exact solution of 
the diffusion equation. 



p(x,0) = 5(x), 



(13) 



the solution on a line can be obtained by Fourier trans- 
form method [4: 



p(x,t) = 



1 



exp 



{x - vtf 
ADt 



(14) 



Similarly, for the circle problem, given the initial condi- 
tion 



P (0,o) = «y(0), 



we obtain 
p((f>,t) = 



, ^ exp 



tot + 2m:) 2 
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In deriving Eq. H16|) . we have used the identity 
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f(x) dx= / f( x + 2n7T ) dx 



(17) 



for a localized function f(x), and we have treated the 
"circle problem" as an "infinite-folded line problem" . 



Furthermore, the "center of mass" now is written 
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exp(— Dt + iu>t) 



(18) 



which is equivalent to Eq.(32) of Ref. [I| if we define 
6 2 2tt 2 2tt 

Dt = —T = —rT, Ujt = 6T= —T. (19) 

2 M 2 M v 1 



Here r and 9 are defined as 
t s 
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Now we compare the results with those in Ref. [l|. 
Fig. d shows the results from Eq. {23), tEl, and J20J) 
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for the cases M = 30 and M = 60. As one can see, they 
indeed reproduce the results of Ref. in the large M 
limit. The parameter N does not appear here because 
the motion of each particle is independent in our model. 

Note that although for a small t the expression l|16|) 
is good enough to be a fast convergent series, however, 
when t becomes large, Eq. (116(1 converges slowly. In this 
situation we use a more accurate expression for p: 
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u>t)], (21) 



n=l 



which can be derived from l|16|) using the Poisson sum- 
mation formula: [f| 



E 



f(na) 



2tt 
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Here /(x) is a localized function, and 
1 



2vr 



/(*) e- 



ri:i; 



(22) 
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is its Fourier transform. 

We now consider some solvable generalizations of 
Eq. Q and JSJ. Note that the ratio between D and 
a; in our model is fixed: 

------ (24) 

to 2 M v 7 

To relax this restriction, we modify our urn model by 
assuming that at each time step the picked ball can have 
probability p to be put into the next urn and probability 
q = 1 — p to be put into the previous urn. Hereafter we 
call this modified model the pq-model. The pq-model is 
also solvable by using methods like those used in Ref. 
. The continuous limit of the pq-model can be derived 
from the recurrence relation for /,;: 



p_ 

N 



f, 
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where 
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to both sides of Eq. (|25() . it becomes 
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Now, define 
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then we get a continuous equation of the form , with- 
out the restriction Ij24(l . One special case is p — q = 1/2, 
which has a zero drifting velocity, and the evolution of 
the system is governed by pure diffusion process — the 
random walk. 

For another generalization we assume that the drifting 
velocity v varies with time, that is 



dp . .dp d 2 p 
_+„(*)_ =2 j_ 

Defining x{t) as the time integral of v{t): 



x(t) 



v{t')dt', 



and adopting the initial condition (|13|) . then 



p(x,t) 



V^wDt 



exp 
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Similarly, for the diffusion equation on a circle with a 
time-dependent uj(t) and initial condition 115(1 : 



(32) 



and the solution is 



p{<t>,t) = 



Here 



, ^ exp 



<t>{t) = / u(t')dt'. 
Jo 



<f>{t) + 2nvr) 



ADt 



(33) 
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The reason for why v and u) can freely vary with time 
relies on Ea. (|28|l . Recall that in our original multiurn 
Ehrenfest model or the pg-model both the time interval 
between two steps and the distance (angle difference) be- 
tween two urns are undefined. Thus in deriving the con- 
tinuous limit of these models we do not have to adopt a 
constant At at each step or a fixed Ax {AO) between two 
neighboring urns. If we relax the restriction in Eq. I|28|) 
and modify them to At s and Ai, (A#j), then the con- 
tinuous limit of these quantities lead to v{t) or ui{i). 

It is interesting to note that the solutions for the dif- 
fusion equation (|10l) can be used to find the wave func- 
tion or Green's function of some time-dependent quan- 
tum mechanical problems |7J . The main idea is to define 
a transformation appropriately between the parameters 
used in the diffusion equation I|1U|) or Q and © and 
those used in the corresponding Schrodinger equations. 
For instance, consider a quantum point particle of charge 



4 



q and mass m moving under the influence of a vector po- 
tential A(t) @]: 



ih 



dtp 
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here we have assumed that A(t) is a function of time t 
only. Rewriting Eq. i|35|) as 
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and executing the transformation 



tp = exp 
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qA(t>) 
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we find 



(36) 



(37) 



(38) 



Comparing Eq. 138|l with (|10|l . we find that they can 
be transformed to each other by the substitution: 
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2m ' 



qA 

) , P 



(39) 



To be more specific, consider the case that the particle 
moving on a circle of radius 1 . Suppose the circle is lying 
on the xy-plane and centered at (x,y) = (0,0). The 
vector potential can be chosen as A(t) = A(t) and is 
generated by a time-dependent magnetic flux $(£) tube 
going through the origin and pointing along the z-axis 



A(t) = A(t) e<t, 



2tt 



Choosing the initially condition as 

V#,o) = V#,o) = <m 



(40) 
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then 
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Here 



U{t) = exp 



2ihm .Iq \ < 



qA{?) 



d£ 



(43) 



and 



Mt) = -— A(t')dt' 
mc J Q 



(44) 



Note that the ip in Eq. H42fl is nothing but the Green's 
function G((p, cpo', t, to) for the quantum particle with 
(j) = to = 0. If A(t) = 0, Eq. g2J> gives the well known 
results |j : 



G{cj>;t) 



m 



exp 



(<p+2nir) 2 
(2iht/m) 



(45) 



and 
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for small and large t, respectively. 

In conclusion, we have derived the continuous limit of 
a mutiurn Ehrcnfcst model, which is a diffusion equation 
with a drifting velocity term. Solving the equation gives 
us the correct time evolution behavior of the ball distri- 
bution. A transformation was introduced, which changes 
the solution of the diffusion equation to the correspond- 
ing solution for the problem of a quantum particle moving 
under the influence of a time-varying magnetic field. 
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